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Abstract 



We propose a general scheme for separation of variables in the integrable Hamilto- 
nian systems on orbits of the loop algebra sl(2, C) x V(\, A^ 1 ). In particular, we illus- 
trate the scheme by application to modified Korteweg — de Vries (MKdV), sin(sinh)- 
Gordon, nonlinear Schrodinger, and Heisenberg magnetic equations. 



Introduction 

y—i ■ 

Let us make a brief review of the problem, 
j — . | After the fundamental paper [23], B. Dubrovin (see [8]) proposed a separation of vari- 

ables for finite gap KdV system. B. Dubrovin shows that the poles of an appropriately 
normalized Baker — Akhiezer function for the auxiliary linear spectral problem are the 
separation variables. The new variables evolve on a hyperelliptic Riemannian surface 1Z 
of genus g. The genus coincides with the number of degrees of freedom of the finite gap 
phase space. 

In the papers [14,18,20,25] a separation of variables is realized for sin-Gordon equation, 
nonlinear Schrodinger equation, and the classic Thirring model. The case of sin-Gordon 
equation appears to be completely similar to the KdV system. However, the cases of 



nonlinear Schrodinger equation and Thirring model have a distinction: the number of 
degrees of freedom is greater by one than the genus of the corresponding spectral curve. 
Here the papers [14, 18, 25] suggest the separation of variables on a reduced phase space. 
Later, the complex Liouville torus of nonlinear Schrodinger equation was proven to be the 
generalized Jacobian of a singular Riemannian surface (see [24]). 

The ideas of the early papers on the integration of finite gap systems were generalized 
by E. Sklyanin [27,28] and partly extended to the quantum integrable models [29,30]. 

At the beginning of the 90s a new technique of separation of variables appeared that 
effectively uses bi-hamiltonian, or multi-hamiltonian, properties of integrable systems, 
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see [3,5-7] and [10-12,22]. The main result in this direction is the diagonalization of 
recursive Nijenhuis operator. In the papers [11,22] the method is applied to KdV and 
Boussinesq hierarchies, and classical finite-dimensional systems. 

The papers [9,31,32] investigate the connection between the problem of separation of 
variables and the parametrization of compact tori by symmetric products of Riemannian 
surfaces. According to [31,32], if a change of variables reduces Liouville 1-form to a sum 
of meromorphic differentials on the corresponding Riemannian surface, then we say that 
the new variables are the separation variables. 

We propose a method of separation of variables for integrable Hamiltonian systems that 
is connected with the orbit structure of affine Lie algebras. The fact that finite gap phase 
space of an integrable soliton hierarchy has an orbital structure was established in [15, 16]. 1 
The Hamiltonian systems in question obey the equations of Lax type, and hence the 
separation variables are points on the corresponding spectral curve. Note that such systems 
are multi-Hamiltonian, which connects our results with the results of [10-12,22]. 

In Sections 1 and 3 we reproduce the key results from [15, 16] about finite gap phase 
spaces for integrable equations as orbits of loop algebra. We illustrate our scheme by 
the examples of modified Korteweg-de Vries (MKdV) system, sin(sinh)-Gordon equation, 
nonlinear Schrodinger equation, and Heisenberg magnetic chain. 

This paper is organized as follows. Sections 1 and 2 are devoted to MKdV system and 
sin(sinh)-Gordon equation. In Section 1 we construct adjoint Poisson spaces and define 
the orbits regarded as phase spaces for MKdV system and sin(sinh)-Gordon equation. 
The construction is discussed in more detail in [4]. In Section 2 we describe the scheme 
for separation of variables and illustrate it by application to MKdV system and sin(sinh)- 
Gordon equation. We show that the separation of variables is achieved on both orbits 
simultaneously. In Sections 3 we construct adjoint Poisson spaces and define the orbits 
regarded as phase spaces for nonlinear Schrodinger equation and Heisenberg magnetic 
chain. In Sections 4 and 5 we similarly consider separation of variables for nonlinear 
Schrodinger equation and Heisenberg magnetic chain, accordingly. 

Acknowledgements. The authors are grateful to participants of the scientific seminar 
'Integrable Hamiltonian Systems and Solitons' T. Skrypnyk, D. Leikin, N. Yorgov for useful 
remarks and discussions. 



1 Phase spaces for MKdV system and sin-Gordon equation 
as orbits in sl(2, C) <g> V(X, A" 1 ) 



First, let us recall some constructions from [15,16]. Take the algebra sl(2,C) with the 
basis 



lr The results of [15, 16] are partially covered by [13]. However the authors of [13] took no notice of 
the remarkable duality between pairs of soliton equations: MKdV and sin-Gordon equations, KdV and 
Liouville equations, nonlinear Schrodinger and Heisenberg magnetic equations, etc. The duality is evident 
if one uses the orbital approach. The pairs of dual equations have common Liouville torus and separation 
variables. Sometimes, there exists a gauge equivalence between the equations of a pair, and the equivalence 
extends to the total infinite phase space [26]. 
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Suppose V(\, A 1 ) is the algebra of Laurent polynomials in A. Denote by g the algebra 
sl(2, C) (g> V(X, A -1 ). Then 

R 2rn = x m H ^ x 2m+l = ynj^ y 2m+l = yrn+ly ^ ^ 

is a basis in g. 

Consider the operator 

d = 2A-^- + &d H ; 
aA 

we call it the operator of principal grading. It is easy to prove that the basis elements (1.1) 
are the eigenvectors of d. We call the eigenvalues of d the degrees. The superscripts in the 
lefthand sides of (1.1) indicate the corresponding principal degrees of the basis elements. 
By di, I G Z, denote an eigenspace of principal degree I. It is evident that 

g 2m = span^tf 2 " 1 }, g 2m+1 = span c {X 2m +\ y 2m+1 }. 

Decompose g into two subalgebras 

9+ = ^2qi, Q- = ^2ei, = 0+ + 0-- 

1^0 l<0 

Further, consider the ad-invariant bilinear forms 

(A(\),B{\)) k = res\- k - 1 TrA(\)B(\), A(X), B(X) eg, keZ. (1.2) 
We use the forms to define the spaces dual to g+ and g_. 
Example 1. Let k = — 1. We have 

(0-)*=5 + +0-i, (fl+r = x>. (i.3) 

where (g_)* and (g+)* contain only the nonzero functionals on g±. 
Example 2. Let fe = N ^ 0. Then 

(fl-)* = Yl 0i ' ®+)* = Yl 0l - 

l^2N+l K2N 

Fix N ^ 0. Consider M ff+1 C g, where an element //(A) G M Ar+1 has the form 

™-($ -^>) 

with 

AT AT+1 AT+1 

a(A) = J] A m a 2m , /3(A) = £ A— 1 ^ 2m _ 1 , 7 (A) = £ A m 72m _L 

m=0 m=0 m=0 

We call M N+l the N-gap sector of g, or shortly the /imie gap sector. 



356 



J Bernatska and P Holod 



Because the factor-algebra Q-/ Yli<-2N-4Si acts effectively on M N+1 , the coadjoint 
action of g_ with respect to the form ( , )_i is well defined on M N+l . The same is true 
for the coadjoint action of g + with respect to the form ( , indeed, the factor-algebra 
S+/I]^2Ar+2S« acts effectively on M N+1 . 

Let C(M N+l ) be the space of smooth functions on M N+1 . For all f u f 2 G C{M N+1 ) 
define the first Lie-Poisson bracket by the formula 

m,n=0 a,b=l " l n 

where 

Cl-^WMC 1 ^' 1 ])-!. 

}J m Z m Y m Z m X m 

H m = a m , fi m = (3 m , fi m = 7 m . 

With the same notation, define the second Lie-Poisson bracket by the formula 

m,n=0 a ,b=l ^ m rn 

where 

P2 n {N) = (MX), [z- m+N ,z- n+N ]) N . 

One can see that the functions P2N+1 and 72N+1 annihilate the bracket (1-4) 

{/W,/}i = 0, {72JV+i,/}i=0 for all / G C(M N+1 ). 

Thus, we can assume without loss of generality that 

P2N+1 = I2N+1 = const (1.6) 

and restrict the bracket (1.4) to the subspace M^+ l C M N+1 with the constraints (1.6), 
clearly, dimM^" 1 = 3(iV + 1). The first Lie-Poisson bracket is nondegenerate on M^ 1 . 
We use the set 72m-i, $i m -\, c«2m, m = 0, 1, . . . , N, as coordinate functions in M^ 1 . 
We call the fixed coordinates P2N+1, 72N+1 the external parameters. 

We see that, on one hand, M^ 1 C ($-)* with respect to ( , see Example 1, and, 
on the other hand, M^ 1 C (5+)* with respect to ( , )jv, see Example 2. 

In addition to the brackets (1.4) and (1.5), one can define N intermediate brackets with 
the Poisson tensors 



PT{k) = </2(A), [Za m+k , Z~ n+k ]) k , k = 0, . . . , N - 1. (1.7) 
Now, consider the ad*-invariant function 

/(A) = - det /Z(A) = h-iX- 1 + h + --- + h 2N +iX 2N+1 . 
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Then we have 

K = ^ («2m«2n + 72m-l#2ra-l) , V = -1, 0, . . . , 27V + I. (1.8) 

The Kostant-Adler scheme [1] implies the following assertions. 

Proposition 1. All functions h u , v = —1, 0, . . . , 2N + 1 determined by (1.8) mutually 
commute with respect to the brackets (1.4), (1.5), and f/ie intermediate brackets with the 
Poisson tensors (1.7). 

Proposition 2. TTie functions h u , v = N, . . . , 2N, are functionally independent and 
annihilate the bracket (1.4). 

Consider the algebraic variety C M^j" 1 defined by the set of equations h v = c„, 
u = N, . . . , 2N, where c v are arbitrary fixed complex numbers. is an orbit of coadjoint 
action of the subalgebra g_ and dimO^ = 2 (TV + 1). 

Proposition 3. The functions h u , v = —1, . . . , N — 1, are functionally independent and 
annihilate the bracket (1.5). 

Consider the algebraic variety O2 C defined by the set of equations h v = c u , 

v = — 1, . . . , JV — 1, where c v are arbitrary fixed complex numbers. is an orbit of 

coadjoint action of the subalgebra g + and dimO^ = 2 (TV + 1). 

It is obvious that the orbits and O2 are the symplectic leaves with respect to the 
first and the second Lie- Poisson brackets, accordingly. 

Further, the functions h-\, ho, /tjv-i, regarded as Hamiltonians with respect to 
the first Lie- Poisson bracket, generate non-trivial flows on 

du a 

-f^H/CMi. ^ = -1,0, ...,7V- 1. (1.9) 

The equations (1.9) can be written with the help of the second Lie-Poisson bracket and 
the functions fajv, . . . , h 2 N regarded as Hamiltonians. Namely (see [16]), one has 

{Vm>h v }l = -{Hm,K+N+l}2- 

Proposition 4. The system (1.9) reduced to the orbit is equivalent to the finite gap 
complex MKdV hierarchy. 

The system (1.9) reduced to the orbit O2 is equivalent to finite gap sin(sinh)- Gordon 
equation. 

Below we give the outline of the proof which may be found in full detail in [4] . 
First, rewrite (1.9) in matrix form 

^ = [V 2 W+i,M(A)] = [£(A),VrU (1-10) 

where 

V1/1 = y ( H- 2m - 2 + y- 2 ™- 1 + x- 2 ™- 1 

^ Q \Oa 2 m OP2m-l Oj2m-l 

y ^ _ ( ®h jj~2m+2N _j_ y-2m+2AT+l _j_ ^h ^-2rn+2N+l\ 

2 ^ \^2m d/3 2m -l dj2m-l J' 
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The hamiltonian flows along r v and t v > commute, which implies the compatibility condition 
in the form of zero curvature equations. In particular, assigning tn~i = x, Tjy-2 = t, we 
obtain 



dV 2 h 2N dV 2 h 2 N-i . rv7 , „ , , n 
at ox 



where 



v 2 h 2N =( x a ™ 



V2/12N- 



«2JV-2 + AQ2AT P2N-I + A/?2Af+l 

,A72at_i + A 2 72at+i -(a 2 Ar_ 2 + Aa 2 iV), 

Recall that 02N+1 an d 72N+1 are the fixed external parameters. 
The reduction of (1.9) onto the orbit gives the equation 

8a2N d&2N-2 



(1.11) 



<9i dx 

equivalent to MKdV equation with respect to the function a2jv(x,t) = u(x,t). Indeed, 
reducing the equation (1.9) as v = N - 1 to the orbit Of we obtain 

I / da 2 N 2 



1 ( , <9a 2 Af 2 



^ = l C27V + - H ' (L12b) 

«2AT-2 = 4 ^ - fa2 ~ 2 °2W + 2 C2Ara2Af^ • (1.12c) 

It is readily seen that combining (1.11) and (1.12c) we get the complex MKdV equation. 
Two real subalgebras su(2) and su(l, 1) = sl(2, R) of sl(2, C) give rise to two real MKdV 
equations (the so-called ±MKdV). 

In the same time, the reduction of (1.9) onto the orbit leads to sin(sinh)-Gordon 
equation. Let O2 U g_i be the base for the orbit O2 '■ The 1-parameter subgroup Go = 
expgo parametrizes the base in a natural way 

7-1 = Vh^~ie u , P-i = ^/h^ie~ u . 

Then the equations (1.9) imply 

a2N= 2{L U > (Li3) 

where x = tjv-i as above; the corresponding flow is called stationary. The Hamiltonian 
hjsj gives rise to an evolutionary flow. In the case of the subalgebra st(2, M) we have 

^ = smhu. (1.14) 

Combining (1.13) and (1.14) we obtain sinh-Gordon equation. 

In the case of the subalgebra su(2) we have to assign 02™ = io2 m , 02™ £ M, and 
72m-i = -#L-i> therefore /3 2 at+i = 727V+1 = ib , 7-1 = -*re m , = -ire" m . Then we 
come to sin-Gordon equation 

d<2u A L 
— — — = 4rosmu. 
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2 Separation of variables for MKdV and sin(sinh)-Gordon 
equation 

Definition 1. Suppose we have the variables (Xk,Wk), k = 1,. . . ,N + 1, such that 

(i) they are quasi-canonically conjugate, that is 

{A fc , wi}i = f{\k)Ski, {A fe , A/}i = {w k , wi}i = 0, 
where /(A) is an arbitrary smooth function; 

(ii) they reduce Liouville 1-form 2 to a sum of meromorphic differentials on the corre- 
sponding Riemannian surface. 

We call (Afc, Wk), k = 1, . . . , N + 1, separation variables. 

Consider the orbit Of , dimO^ = 2(N + 1). One can parameterize the orbit using 
any subset of 2(N + 1) variables from {«2m, fiim-i, 72m-i}j m = 0,1, . . . ,N. The most 
natural way to obtain the parameterization is to eliminate one of the subsets {ftm-i} or 
{l2m~i}- The reason for this is the nilpotency of the basis elements that correspond to 
the subsets. 

Note that the correspondence between the elimination variables, which we chose to 
parameterize the orbit, and the nilpotent elements of the basis of the algebra is a crucial 
feature of our scheme and applies to all examples. 

We chose to parameterize the orbit Of by the variables {^m-i, ct2m}, m = 0, 1, . . . , N, 
that is we eliminate the set {fom-i}- From the orbit equations we find 

N+l 

fom-l = ^2 ( r+ )mj( c A r +J ~ A N+j), m = 0, ... N+l, C 2 N+1 = P2N+H2N+I, (2.1) 
j=0 



where 



72JV+1 12N-1 
72Af+l 








7i 

73 



7-1 
7i 



72JV+1 72AT-1 
72N+1 



and 



A,. 



E 



m+n=v, 



Now, using the parameterization (2.1), we find expressions for the Hamiltonians h-i,ho, 
■ ■ ■, hisr-i 

N+l 

K-i = r- m (r + )-;.( C7V+ , - A N+3 ) + A n _ x , n = 0,...N, (2.2) 



where 



m,j=0 

7-1 
71 





7-1 



72JV-1 72iV-3 



7-1 



2 We call Liouville 1-form if dQ — u, where uj is a symplectic 2-form. 
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Note that the expressions (2.2) are linear in c v , v = N, . . . , 2N + 1. 

Clearly, one can obtain an analogous parametrization of the orbit by the set 

{«2m, 72m-l}> m = 0,...,N. 

To proceed we need to define the characteristic polynomial 

Q(x, A) = det(/x(A) - x-i"), 

where / denotes 2x2 identity matrix. By the substitution x = tuA -1 the equation 
Q(x, A) = becomes transformed into the standard equation of a hyperelliptic curve of 
genus N + 1 

P(w/, A) = \ 2 Q(w\-\ A) = w 2 - A(/i_i + /i A + • • • + /i 2 jv+iA 27V+2 ) = 0. (2.3) 

Recall that on the orbit 0^ we have h„ = c v , v = N, . . . , 2N. Denote by (wk, A&) a 
root of P(w, A) on the orbit, that is 

w 2 = A fc (/i_i + Mfc + • • • Ziiv-iAf + Cjv Af +1 + Cjv+1 Af + 2 + • • • + c 2N+1 X 2 k N +2 ). (2.4) 

We proceed to show that the set {(u>fc, Ajt)}, A; = 1,...,JV+1, defines another parametriza- 
tion of the orbit . We have to find the explicit relation between the sets {(u>i, Ai), . . . , 
(wn+i, Aat+i)} and {a , a 2 , . . . , a 2N , 7-1,71, . . . ,727V-i }• 

Solving (2.4) for the Hamiltonians h-i, ho, ... , h^-i one gets 

fc-i = ^[Wi - ^(A^ 1 ) c 2 ^ + i^i(A 2Ar+2 )] 

to = ^2 (£) - ^(A^ 1 ) c 2N+1 W 2 {\ 2N + 2 )] (2 5) 

/ijv-i = ^[Wjv+i (^) - c JV W) v+ i(A JV+1 ) c 2N+1 W N+1 (\ 2N+2 )}. 

where W = Yi(\ ~ i s Vandermonde determinant of Ai, X 2 , . . . , Ajv+i- By Wi(f(X,w)) 
we denote the determinant of Vandermonde matrix with the i-th column replaced by 
{f(M,w 1 ),...,f(X N+1 ,w N+1 )) • 

On the orbit the formulas (2.2) and (2.5) define the same set of functions. We see that 
both (2.2) and (2.5) are linear in c v , v = N, . . . , 2N + 1. As {c u } is the set of independent 
parameters one can equate the corresponding terms. Namely, we get 

^l,-""f I m = 0, . . . JV. 

727V+1 W 

This implies that the set {Afc} is, in fact, the set of roots of the polynomial 7(A) 
7(A fe ) = 0, 

while the variables {wt} satisfy the equalities 

w 2 = A 2 (a 2 (A fc ) - 7(A fc )0(A fc )) = A 2 a 2 (A fc ), k = 1, . . . , N + 1. 
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Theorem 1. Suppose the orbit C?i has the coordinates (a2 m , 72m-i)> m = 0, 1, . . . , N, 
as above. Then the new coordinates (Afc, Wk), k = 1, . . . , N + 1, defined by the formulas 

7(A fc ) = 0, w k = e\ k a(X k ), where e 2 = 1, (2.6) 

have the following properties: 

(1) a pair (wk, Afc) is a root of the characteristic polynomial (2.3). 

(2) a pair (Xk,w k ) is quasi- canonically conjugate with respect to the first Lie-Poisson 
bracket (1.4): 

{Afc, Xi}i = 0, {X k , wi}i = eX k S k i, {w k , w{\i = 0; (2.7) 

(3) the corresponding Liouville 1-form is 

= y^eX^Wk dX k - 



Proof. (1) The assertion is a direct consequence of (2.3) and (2.6). 
(2) It is evident that 

{A fe ,A z }! = 0. 

Indeed, since Afc, k = 1, . . . , N + 1, depend only on 72m-ii = 0, . . . , N, and 72m- 1 
mutually commute, Afc also mutually commute. 
Let us calculate the bracket of Afc and wi 

s-^ ( dX k dwi dX k dwi \ 

{Xk,Wi} 1 = y |72m-l,Q!2n}l- 

^ V^72m-i oa 2n oa 2n oj 2m -i ) 
From (2.6) we have 

-^- = 0, &Xk = A " — =cA n+1 . (2.8) 

Further {7 2m _i,a 2 „}i = -7 2 (m+n)+i wh en m + n < N and {72m-i, «2n}i = when 
m + n ^ N. Thus, we obtain 

E £^A[ l+1 72(m+n)+l . . , , 

rx -, m+n<N eXi 7(A«J-7(Afc) 

7 (Afc) 7 '(Afc) A/ - A fc 

As fc / / it is evident that {X k ,wi}\ = while 7(A/) = j(Xk) = 0- As /c = / we get 

/x , ,. eAj 7(^)-7(A fc ) 
{A fc ,wfc}i= bm = eA fc . 

Ai^A fe 7'(A fc J A; - Afc 

Thus, 

{Afc, = eAfc^. 
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Let us calculate the bracket of w k and wi 



r i ST- f dw k dwi dw k dwi \ 

{w k ,Wi}i = > ^ ^ ^ |72m-l,a2n}l- 

^ V^72m-i oa 2n oa 2n oj 2m -i J 
From (2.6) follows that 



r , , , <9A fc 



e [a(Afc) + AfcQ'(A fe )] 



_ i y * tv J 1 ' *fi< V ft / I O 

»72m-l <?72ni-l 

then, using (2.8), we obtain 

r , fA,[a(A fc ) + A fc o'(A fc )] A fc [a(Aj) + Aja^Aj)]^ 7(Aj) — 7(A fe ) 

{W k , Wl\\_ — 

hence 



y(A fc ) y(A,) y A z -A fc ' 

= 0. 



(3) From (2.7) it follows that Liouville 1-form on the orbit is 

VL-i = y^ eX^Wk dX k . 
k 

The reduction to Liouville torus is done by fixing the values of Hamiltonians h—i, 
ho, /ijv— i- On the torus w k is the algebraic function of A& due to(2.3). After the 
reduction the form 0_i becomes a sum of meromorphic differentials on the Riemann 
surface P(w, A) = 0. ■ 

The next theorem is proven similarly. 

Theorem 2. Suppose the orbit 0^ has the coordinates (a2 m , 72m-i)> m = 0,1,..., N. 
Then the new coordinates (Xk,Wk), k = 1, . . . , N + 1, defined by the formulas 

7(A fc ) = 0, w k = eX k a(X k ), where e 2 = 1, 

have the following properties: 

(1) a pair (wk, X k ) is a root of the characteristic polynomial (2.3); 

(2) a pair (Xk,Wk) is quasi- canonically conjugate with respect to the second Lie-Poisson 
bracket (1.5): 

{X k ,Xih = 0, {X k ,w l } 2 = -eX k v+2 5 k i, {w k , Wl } 2 = 0; (2.9) 

(3) the corresponding Liouville 1-form is 

tt N = - ^ £ Afe (JV+2 W dA fc . 
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Let us summarize our scheme of obtaining the separation variables. First, we param- 
eterize the orbit by eliminating a subset of group coordinates corresponding to nilpotent 
basis elements. Next, we restrict the curve P(w, A) = onto the orbit, where P(w, A) is 
the characteristic polynomial 

P(w,\) = det (fi(X) - w I), 

I is identity matrix. We use the set {Afc, Wk}, k = 1, . . . , N + 1, where P{wk, A&) = 0, to 
define another parametrization of the orbit. Then, we equate expressions for Hamiltonians 
in the coordinates of two parameterizations of the orbit in order to obtain the link between 
the two sets of orbit coordinates. Finally, the set {A&, w^} is the set of separation variables. 

Further, in Sections 4 and 5 we apply the scheme to nonlinear Schrodinger equation 
and Heisenberg magnetic chain. 

3 Phase spaces for nonlinear Schrodinger equation and 
Heisenberg magnetic chain as orbits in s((2, C) <g> V(z, z~ x ) 

Here we use the construction from Section 1 with homogeneous grading. That is, 

X 1 = z l X, Y l = z l Y, H l = z l H (3.1) 

be the basis in g ~ sl(2, C) (g) V(z, z' 1 ). 

Note the well-known fact that the Lie algebra from Sections 1-2 can be realized as the 
subalgebra of sl(2, C) ®T>{z,z~ l ) invariant with respect to an automorphism of order 2, 
see [17], [19], [21]. 

By gi, I G Z, denote an eigenspace of homogeneous degree I. It is evident that 

Ql = span c {X', Y\ H 1 }. 
Decompose g into two subalgebras 

9+ = ^2qi, 8- = ^2m, = ?++?-• 

1^0 KO 

Use the same ad-invariant bilinear forms (1.2) to define the spaces dual to g+ and g_. 
Fix N ^ 0. Consider M N+1 C g, where an element ju(z) G M N+1 has the form 

\-y(z) -a{z)J 

with 

N+l N+l N+l 

m=0 m=0 m=0 

As above, we call M N+l the N-gap sector of g, or shortly the finite gap sector. 
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For all /i,/2 € C(M N+l ) define two Lie-Poisson brackets 

, d_h df 2 



N+1 3 At At 



and 



m,n=0 a, 6=1 



{/i,M2= E E p r^+i)^^, (3.3) 



m,n=0a, b =l d ^ 9 ^ 



where 



PTi-i) = {Kz)AZa m -\z b n - 1 ])-u 

PT{N + 1) = <£(*), [^- m+Ar+1 , zr n+iV+1 ])Ar + i, 

Mm = Mm = Ati> Mm = 7m- 

One can see that M N+1 C ($-)* with respect to ( , )_i and, in the same time, 
M N+1 C with respect to ( , )jv+i- 

Next, introduce the ad*-invariant function 

I(z) = - det Ji(z) = h + h 1 z + --- + h 2N +2Z 2N+2 , 

where 

K= E («™«" + 7m/3n) , v = 0, 1, . . . , 2N + 2. (3.4) 
m+n=;/ 

One can easily prove that the functions ajv+i, Pn+i, Jn+i annihilate the bracket (3.2). 
In order to obtain nonlinear Schrodinger equation we have to assign 

(3 N+1 = 7at + i = 0, a N+ i = const / 0. (3.5) 

After the restriction of the bracket (3.2) to the subspace C M N+l with the con- 

strains (3.5) we get dimM^ 1 = 3(iV + l). The bracket (3.2) is nondegenerate on M^ 1 . 
We use the set j m , (3 m , a m , fn = 0, 1, . . . N as coordinate functions in M^ 1 . We call the 
fixed coordinate ajv+i the external parameter. 

On the other hand, the functions ajv+i, @n+i, 1n+i commute with all Hamiltonians h u , 
v = N + 2, N + 3, . . . , 2N + 2, with respect to the bracket (3.3) and give rise to nontrivial 
flows on M N+1 , therefore are Hamiltonians. That is, the bracket (3.3) is considered on 
M N+1 , dimM Ar+1 = 3(N + 2), and the set -y m , /3 m , a m , m = 0,1, ... N + 1 serve as 
coordinate functions in M N+1 . 

The following assertions are immediately derived from the Kostant-Adler scheme [1]. 

Proposition 5. All functions h v , u = 0, 1, ... , 2N + 2 determined by (3.4) mutually 
commute with respect to the brackets (3.2) and (3.3). 
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Proposition 6. The functions h u , u = N + 1, . . . , 2N + 1 are functionally independent 
and annihilate the bracket (3.2). 

Consider the algebraic variety C defined by the set of equation h v = c v , 

v = N + 1, . . . , 2N + 1, where c u are arbitrary fixed complex numbers. Oi is an orbit of 
coadjoint action of the subalgebra g_ and dimO^ = 2(iV + 1). 

Proposition 7. The functions h u , v = 0, . . . , N + 1 are functionally independent and 
annihilate the bracket (3.3). 

Consider the algebraic variety 0^ +1 C M N+l defined by the set of equation h u = c v , 
v = 0, . . . ,7V + 1, where c v are arbitrary fixed complex numbers. 0^ +1 is an orbit of 
coadjoint action of the subalgebra g + and dimO^ +1 = 2(iV + 2). 

Further, the functions ho, hi, . . . , hjy regarded as Hamiltonians with respect to the 
bracket (3.2) give rise to nontrivial flows on M^ 1 

du a 

-jPM/CMi. u = 0,l,...N. (3.6) 

The functions hw+2, • • • , h 2 N+2 regarded as Hamiltonians with respect to the bracket (3.3) 
give rise to nontrivial flows on M N+l 

da a 

-f^ = -{/C,^+iV+2}2, v = 0,1, . . . N. (3.7) 

Proposition 8. The system (3.6) reduced to the orbit is equivalent to finite gap 
nonlinear Schrddinger equation. 

The system (3.7) reduced to the orbit 0^ +1 is equivalent to finite gap Heisenberg mag- 
netic chain. 

Here we give the outline of the proof. 

Consider the orbit 0± . Rewrite (3.6) in matrix form. In particular, assigning tn = x, 
tn-i = t, we obtain 

= $(z), Vi/ijv] = [V 2 /i 2 JV+i, Hz)], (3.8a) 

^ = $(z),\7ih N -i] = [V 2 h 2N ,ft(z)}, (3.8b) 

where 

y h _( Z012N+1 + OLN Pn 

2 2N+1 V 7Af -(za 2 N+i + a N )J ' 

/ z 2 a 2 N+i + za N + a N -i z(3 N + /3jv-i 

V z77v + 7Ar_i -{z z a N+ i+ za N + a N -i 

We use the real subalgebra su(2) of sl(2, C), that is assign a m = ia m , 7 m = =F^, then 
the compatibility condition for (3.8) gives 

df3 N d 2 (3 N 2 

2zaAr + i-^^ = —~q^2 *PN\Pn\ — APN">2N, 
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which coincides with nonlinear Schrodinger equation with respect to the function /3n(x, t) = 
if)(x,t) as ci 7v +1 = ^, = a,N = 0, 

Consider the orbit Of +1 . Note, that dimO^ +1 = 2{N + 2) and the set h N+2 , h N+3 , 
h>2N+2 is insufficient to provide Liouville integrability. Since the functions ajv+i, 
/Sjv+i) and 7n+i are in involution with the set /iat+2, ■ • • , ^-2iV+2 with respect to the 

bracket (3.3) one can take any of them as an extra Hamiltonian. Here we chose ccjv+i- 

Rewrite (3.7) in matrix form. In particular, by assigning t n+2 = x and tn+3 = t we 
obtain 

- [V 2 /ijv+2,?(2)] = [£(*), Vi/io], (3.9a) 



dx 
dj2(z 



dt 

where 

-i /" "o A) 



[V 2 /ijv+3,£(*)] = [^),Vi/ii], (3.9b) 



7o -«o 



\7i -«iy \7i -«i 

Let us replace the coordinates a m , /3 m , j m , m = 0, 1, . . . , N + 1, according to the 
formulas 

Gto = i^rni fim = f^m ~ Vmi 7m = — AVra — Vm 1 

Then 

— £ ijk^ m +n-N-l- 

Introduce the vector notation /x m = (//^, // m )*, m = 0, 1, . . . , N + 1. Now the orbit 
is determined by the equations 

Quo, Mo) = -co, 
2(/i ,/Xi) = -ci, 

I] (/*m)/*n) = — CJV+l, 
m+n=Af+l 

where (-, •) denotes the dot product. The equations (3.9) are written in the form 
dfx n 



dx 

dll m 



2[/i ,/*m+i], (3.10a) 
2[/xi,/x m+ i] + 2[/* ,/*m+2], (3.10b) 



0t 

where [•, •] denotes the cross product. 



On Separation of Variables for Integrable Equations of Soliton Type 



367 



By reduction of (3.10) to the orbit 0^ +1 we obtain 



Mi 



1 

2c~o 



Mo, 



duo 

dx 



ci 

+ 7T- MO- 

2c 



Taking into account the compatibility condition 
duo dpi 



dt 



dx ' 



we get 



dt 



1 

2^o 



Mo, 



d 2 Mo 

dx 2 



+ 



ci d^o 

2co dx 



(3.11) 



When ci = 0, the equation (3.11) becomes the well-known classic Heisenberg magnetic 
equation, also called isotropic Landau- Livshits equation. 



4 Separation of variables for nonlinear Schrodinger equa- 
tion 

Consider the orbit , dim Of = 2(N + 1). The most natural way to parameterize the 
orbit, which we already noted in Section 2, is to eliminate the subset {f3 m } (or {7m}), 
m = 0, 1, . . . , N + 1. Then, roots of 7(A) (or /3(A)) give a half of separation variables. 
This way is applied in [18]. However, in the case of nonlinear Schrodinger equation the 
finite gap phase space is determined by the constrains (3.5). That is, the polynomials /3(A) 
and 7(A) have the order N, therefore the set of roots is insufficient to parameterize the 
orbit Of. 

In order to solve this problem we change coordinates. We use coordinates as in [24]. 
Let 

7). ij). *-(Jj Jj). («) 

be the basis in sl(2, C). It is easily shown that [T, S] = S, [T, R] = -R, [S, R] = 2T. Then 
is a basis in g ~ sl(2, C) <8> V(z, z^ 1 ). An element € M^ 1 has the form 



where 



\[r{z) + s{z)\ Hr( z ) - s(z) - 2t(z)] 
] [s(z)-r(z)-2t(z)} -\[r{z) + s(z)) 



N+l N+l N+l 

t(z) = r(z) = z m rm, s(z) = £ z 

m=0 m=0 m=0 



m 

*>r. 
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Note that t m , r m , s m , m = 0, 1, . . . , N + 1 are defined by the formulas: 
-m-i ) _ i = (j2( z ),T- m+N+1 

-m—1 



t m = (j2(z), T _ " t_1 )_i = T- - ■ - ■ - )N+1 , 

r m = @(z), R- m - 1 )^ = ($(z), R~ m+N+1 ) N+1 , (4.2) 
s m = @(z), S-" 1 " 1 )-! = 5- m+JV+1 )AT +1 , m = 0, 1, . . . , N + 1. 

With the new coordinates M^ 1 is defined by constrains 

tN+l = 0, rjy + i = Sn+1 = \] C2N+2 / 0. 

One can see that the subsets {r m } and {s m }, m = 0, l,...,iV + l have nilpotent 
corresponding basis elements. For the reason we chose one of these subsets, namely {r m } 
here, to parameterize the orbit . From the orbit equations we find 

N+l 

r m =^2 ( S+ )mj( c j+N+i - B j+N+1 ), m = 0, . . . , N + 1, c 2 n+2 = r N+1 s N+1 , (4.3) 

j=0 



where 



S + 



SN+1 sn 
s N+ i 













si 

S2 



so 

Si 



sn+1 sn 
s N+ i 



and 



B v — y t m t n . 



m+n=u, 
0^m,n^N 



Using the parameterization (4.3), we find the expressions for the Hamiltonians ho, hi, . . . , /tjy 



N+l 



K = Yl S nm(S + ) m ){Cj+N+l ~ B j+N+l ) + B n , XI = 0, . . . N, 

mj=0 



(4.4) 



where 



S~ = 



so 

Sl 





«0 








_s N Sat_i ... s 0_ 
To proceed we define the characteristic polynomial 
P(w, z) = det(/i(z) —w-l). 



(4.5) 



The equation P(w, z) = has a form of the standard equation of a hyperelliptic curve of 
genus N + 1 



P(w, z) = w 2 - (h + hiz + • • • + h2N+2Z 2N+2 ) = 0. 



(4.6) 



On the orbit we have h v = c v , v = N + 1, . . . , 2N + 1. Denote by (wk, zjS) a root 
of P(w, z) on the orbit, that is 



wl = ho + h x z k -\ h N Zk + Cjv+i4 v+i H c 2 jv+2Z 



JV+l 



(4.7) 
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We proceed to show that the set {(w k ,z k )}, k = 0, 1, . . . , TV + 1 defines another pa- 
rameterization of the orbit . We have to find the explicit relation between the sets 
{(wi, zi), {w N+ i, z N+1 )} and {to, h,..., t N , s ,si,..., s N }. 
Solving (4.7) for Hamiltonians ho, hi, . . . , one gets 

h = hw x {w 2 ) - c N+1 Wi(z N+1 ) c 2N+2 Wi(z 2N+2 )} 

hi = — \W 2 {w 2 ) - c N+ iW 2 (z N+1 ) c 2N+2 W 2 (z 2N+2 )] (4g) 

h N = t\[W n+ i {w 2 ) - c N+1 W N+1 (z N+1 ) c 2N+2 W N+ i (z 2N+2 )}, 

where W and Wi(f(z,w)) denote the same as in (2.5). 

On the orbit the formulas (4.4) and (4.8) define the same set of functions. We see 
that both (4.4) and (4.8) are linear in c u , v = N + 1, . . . , 2N + 2. As {c u } is the set of 
independent parameters one can equate the corresponding terms. Namely, we obtain 

Sm W m+1 (Z N+1 ) 

= 777 , m = 0,l,...,JV. 

SN+1 W 

This implies that the set {zk} is the set of roots of the polynomial s(z) 
s{z k ) = 0, 

while the variables {wk} satisfy the equalities 

w 2 k = t2 (zk) + s(z k )r(z k ) = t 2 (z k ), k = 1, . . . , N + 1. 

Theorem 3. Suppose the orbit has the coordinates (t m ,s m ), m = 0,1, . . . , N , as 
above. Then the new coordinates {z k ,w k ), k = 1, . . . ,N + 1, defined by the formulas 

s(zk) = 0, Wk = et(zk), where e 2 = 1, (4-9) 

have the following properties: 

(1) a pair (w k ,z k ) is a root of the characteristic polynomial (4.6). 

(2) a pair (z k ,w k ) is canonically conjugate with respect to the Lie-Poisson bracket (3.2): 

{z k ,zi}i = 0, {z k ,wi}i = e5 k i, {w k , wi}i = 0; (4.10) 

(3) the corresponding Liouville 1-form is 

= ^eujfe dz k . 
k 

Proof. (1) The assertion is a direct consequence of (4.6) and (4.9). 
(2) It is evident that 



{z k ,z t }i = 0, 
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since Zk, k = 1, . . . ,N + 1 depend only on s m , m = 0, 1, . . . , TV and s m mutually commute. 
Let us calculate the brackets {zk,wi}\ and {wk,wi}\. From (4.9) we have 

dz k _ dz k _ z% dwi _ n dwi _ , dz t 



dt n os m s'(Zk) ot n as m ds m 

Further {s m , t n }i = —s m+n when m + n ^ N and {s m , t n }\ = when m + n > N . Thus, 
we obtain 

r _ g s(z k ) - s(zi) _ ( t'(zk) _ t'(zi) \ s(z k ) - s(zi) 

fc ' ' 1 s'(zk) Zk-zi ' k, l i y s '( Zk ) s'(zi)J z k -zi 

Thus, 

{z k ,wi}i = eSki, {w k , wi}i = 0. 
(3) From (4.10) it follows that Liouville 1-form on the orbit Of is 

Vl-i = y^gwfc dz k . 
k 

The reduction to Liouville torus is done by fixing the values of Hamiltonians ho, hi, ... , 
h 7v . On the torus Wk is the algebraic function of Zk due to (4.6). After the reduction the 
form 9-i becomes a sum of meromorphic differentials on the Riemann surface P(w, z) = 0. 



Given a set of pairs (z k ,Wk), k = 1, . . . , N + 1, one can find the set (t m ,s m ), m = 
0, . . . , N, such that the equations (4.9) are satisfied. Thus, we can define a homomorphism 

c 2N+2 ^ Q N (4n) 

that takes the set of pairs (zk,Wk), k = 1, . . . , N + 1, to a point of the orbit Of . Since 
the orbit Of is topologically trivial, this implies that the map (4.11) is global. 

After the reduction (4.11) turns into the map that takes the (N+ l)th symmetric power 
of Riemann surface to Loiuville torus: 

Sym{7e x K x • • • x 11} ^ T N+1 . 



5 Separation of variables for Heisenberg magnetic chain 

Consider the orbit Of +2 , dimO^ +2 = 2(iV + 2). We chose to parameterize the orbit 
Of +2 by the variables {^ m , a m }, m = 0, 1, . . . , N + 1, that is we eliminate the set {/3 m }, 
which corresponding basis elements are nilpotent. 
From the orbit equation we find 

N+l 

An=£ (f-)-J( Cj -Aj), m = 0,...N + l, (5.1) 

3=0 
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where 



7o 
7i 





7o 



7AT 77V- 1 

.7jv+i In 



0' 



70 

71 7o. 



and 



E 



m+n=v, 
0^m,n^N+l 



Now, using the parameterization (5.1), we find expressions for the Hamiltonians h^ +2 , 
hN+3, ■ ■ ■ , h 2 N+2 



N+l 



K+N+l = r «m( r )m)( C 3 " A i) + A n+N+1, Tl = 1, . . . N + 1, 

mj=0 



(5.2) 



where 



7Ar+i 








72 7i 

73 72 

77v+i 



Note that the expressions (5.2) are linear in c„, v = N, . . . , 2N + 1. 
To proceed we use the same characteristic polynomial (4.5) 

On the orbit +1 we have h u , v = 0, 1, . . . , N+l. Denote by (wk, Zk) a root of P(w, z) 
on the orbit, that is 



w'i = c + c\z k H CAr + izf +1 + h N+2 z£ +z H h 2N+2 z z k 



N+2 



y 2N+2 



(5.3) 



The set {{w k , z/.)}, k = 0, 1, . . . , TV + 1 is insufficient to parameterize Let us fix 

ayv+i and 77V+1 regarded as Hamiltonians and consider below the reduced orbit 2 rtd- ^ 
we find an explicit relation between the sets {(u>i, z±), . . . , (wjv+i, zjv+i)} and {ao, ai, . • • , 
ajv, 70, 71, . . . , 7at} we show that the set {(wk, z k )}, k = 0, 1, . . . , N + 1 defines another 
parameterization of the orbit O^j^ 1 . 

Theorem 4. Suppose the orbit has the coordinates (a m ,^ m ), m = 0, 1, . . . , N, as 

above. Then the new coordinates (zk,Wk), k = 1, . . . , N + 1, defined by the formulas 

l(zk) = 0, w k = ea(zk), where e 2 = 1, (5.4) 

/jaue i/te following properties: 

(1) a pair (w^, z^j is a root of the characteristic polynomial (4.6). 

(2) a pair (z^^Wk) is quasi- canonically conjugate with respect to the Lie-Poisson bracket 
(3.3); 



{Zk,Zl}^ 



0. 



{zk,wih 



-ez 



N+2 



{w k ,wi} 2 = 0; 



(5.5) 
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(3) the corresponding Liouville 1-form is 

^N+i = - ^ez^ N+2 ^w k dz k . 
k 

Proof. (1) The assertion is a direct consequence of (4.6) and (5.4). 
(2) It is evident that 

{zk,zi} 2 = 0, 

since z k , k = 1, . . . , N + 1 depend only on 7 m , m = 0, 1, . . . , N and 7 m mutually commute. 
Let us calculate the brackets {z k ,wi} 2 and {w k ,wi} 2 . From (5.4) we have 

dz k dz k 4™ dwi n dwi , dz t 

0, -= — = — 77^, -= — =ez?, t. — =ea(zi)- 



da n ' d"/ m i{z k y da n "' dj m d^ m 

Further {j m , a n } 2 = -j m+n -N-i when m + n < N + 1 and {7™, a n } 2 = when m + n > 
N + 1. Thus, we obtain 

{z k ,wi} 2 - 



{w k ,wi} 2 
whence 



l'{z k ) z k - zi 

1 1 \z»+*>y(z l ) - z^lizk) 



i{z k ) Y(zi)J Z k -Z t 
{z k ,wi} 2 = -z^ +2 5 kU {w k , w{\ 2 = 0. 



(3) From (5.4) it follows that Liouville 1-form on the orbit O^r^ 1 is 
Oat+i = - y^ez^ N+2) w k dz k . 

k 

The reduction to Liouville torus is done by fixing the values of Hamiltonians ho, h±, . . . , 
hjy- On the torus w k is the algebraic function of z k due to (4.6). After the reduction the 
form fi_i becomes a sum of meromorphic differentials on the Riemann surface P(w, z) = 0. 

■ 

Remark 1. When 7^+1 = 0, a N+ i / one can replace a m ,(3 m ,j m by t m ,s m ,r m , m = 
0, 1, . . . , N, according to (4.2). In this case the new coordinates (z k , w k ), k = 1, . . . , N + 1, 
defined by the formulas 

s{z k ) = 0, w k = et(z k ), where e 2 = 1. (5.6) 

and have the same properties as in Theorem 4. 

The results of Theorem 3 and 4 can be summarized as follows. Liouville tori for the 
nonlinear Schrodinger equation and Heisenberg magnetic chain have the same number of 
parameterizing variables z k and each variable belongs to the hyperelliptic curve (4.6) of 
genus g = N + 1. In other words, the common Liouville torus is the Jacobi variety of the 
curve (4.6). 
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